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Abstract 

A complete perturbative expansion for the Hamiltonian describing 
the motion of a quantomechanical system constrained to move on an 
arbitrary submanifold of its configuration space R n is obtained. 



1 Introduction 



The quantomechanical description of constrained systems is extremely im- 
portant in physics and since the early days of quantum mechanics several 
techniques have been developed to deal with this matter. A fundamental 
contribution has been given by Dirac [0]. His idea, geometrical in nature, 
consists in removing the redundant degrees of freedom by the construction 
of a consistent hamiltonian formalism for the constrained classical theory and 
proceeding then to its quantization. Other noteworthy approach have been 
developed by Schwinger, Peierls by using variational arguments and by De 
Witt, Faddeev, Popov by means of lagrangian formalism and of path integral 
techniques. The common feature to all this methods, to which we will refer 
as formal methods, is that the reduction of the dynamics to the constraint 
surface is performed before the quantization of the system. Although in 
many cases this is the only way to proceed, it often introduces non physical 
ambiguities causing some pure quantomechanical effects to be ignored. 

Let us consider for example a particle constrained to move on an arbi- 
trary surface £ embedded in the three-dimensional euclidean space R 3 . The 
reduction of the classical theory is straightforward. Introducing x x ,x 2 co- 
ordinates parametrizing the surface and denoting by g^ ui fi, v = 1, 2, the 
metric induced on E from R 3 , the system is described by the lagrangian 
C = \g^, u x >J 'x 1 ' . Defining the generalized momenta p^ = dCjdx^ we obtain 
the Hamiltonian TL = \g^ u p^p Ul 9^ u denoting the inverse of the metric. The 
quantization of 7i contains ordering ambiguities which are not completely 
removed by the required covariance of the theory. As first observed by De 
Witt ||, in constructing the Hamiltonian operator we are free to add to mi- 
nus one half the laplacian A a term proportional to the scalar curvature R 
of the surface, 

H = ~A + aR. (1) 

Different quantization schemes produces different values for the constant a. 
a may not be unambiguously determined, depending essentially on an or- 
dering choice. On the other hand the reduction of the motion of a particle 
to a surface is by no means an academic problem. Devices producing the 
confinement of electrons on a plane are widely studied in physics, for exam- 
ple in the Quantum Hall Effect's context ||, and we may think to use the 
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same techniques to constraint a particle on an arbitrary surface S. How to 
determine then the constant a? In addressing the solution of the problem it 
is convenient to give up the formal treatment of the constraint, thinking in- 
stead at the physical mechanism producing the confinement of the particle to 
the surface. The analysis of devices used in the Quantum Hall Effect stuff, 
suggests the confinement to be produced by a potential presenting a deep 
minimum in correspondence of the constraint surface and depending only on 
the coordinate normal to it |3], |4], [5]. Moving along this line H. Jensen and H. 
Koppe II have first given a realistic description of the motion of a particle 
on a surface embedded in R 3 . In accordance with Heisenberg's principle the 
confinement causes the particle to fluctuate very strongly in the direction 
normal to the surface so that the spectrum of the system is described in first 
approximation by that of the confining potential. In correspondence of each 
level the effective Hamiltonian describing the motion along the surface may 
then be unambiguously obtained as (see Ref. |J and the discussion below) 

H = -\ A + \ R -\ e > (2) 

where £ is the extrinsic mean curvature of the surface E. The analysis have to 
be completed by considering the interactions between the degrees of freedom 
normal to the surface and those parallel to it. The case of a wire embedded 
in R 3 has also been considered by many authors 

From this simple example we learn that the operation of reducing the 
dynamics and quantizing a constrained system do not in general commute 
and that performing the former before the latter may produce the appearance 
of unphysical ambiguities and the neglecting of contributions connected to 
the extrinsic geometrical properties of the constraint surface. 

In this paper we present a complete perturbative description of a system 
constrained to move on a submanifold of its configuration space R n by a 
confining potential Vc- Generalizing the example of the surface embedded 
in R 3 we require Vc to satisfy two very general conditions 

CI) Vc presents a deep minimum in correspondence of the constraint sur- 
face, 

C2) Vc depends only on coordinates normal to the constraint surface. 

In our scheme the reduction of the dynamics to the constraint surface is 
performed after the quantization of the system. Adapting coordinates to the 



3 



constraint in section 2, a complete perturbative expansion for the Hamilto- 
nian describing the motion of the system is obtained in section 3. In accor- 
dance with Heisenberg principle the zero order term of the expansion takes 
into account the fluctuations of the system in the directions normal to the 
constraint surface. The first order term, already discussed in Refs.f|, [|, §], 
describes the effective constrained dynamics while the rest of the perturbative 
expansion describes the interactions between normal and effective degrees of 
freedom. We want to point out that the constrained quantum dynamics is 
characterized by the whole expansion, the effective dynamics on the con- 
straint surface representing only the leading term. Contrary to what hap- 
pens in the classical description, the explicit form of the potential realizing 
the constraint leaves traces in the effective dynamics and in the spectrum of 
the system and therefore may not be neglected. This is illustrated by two 
example in sections 5 and 6. Section 7 contains our conclusions. 



2 Geometrical Preliminary 

In the sequel we identify the constraint surface with a smooth m-dimensional 
submanifold M of the configuration space R n . Denoting by $ : M — > R n the 
embedding of M in R n and by n 1 (x), n 2 (x), n n ~ m (x), a smooth assignment 
of (n—m) orthonormal vectors normal to M in every point x G M, an adapted 
coordinates frame may be introduced by using coordinates a^,/x = 1, ...,m, 
on M, plus the distances y\i = l,...,n — m, along the geodetics leaving 
M with speed n\ In a "suflicently close" neighbourhood of M the frame 
{x^, y l ; fi = 1, m, i = 1, n — m} is well defined and its relation with the 
Cartesian coordinates r = (r 1 , ...,r n ) of R n is given by 

r = +y i n i {x^). (3) 

It is important to realize that the embedding of M in R n is completely 
characterized by the assignment of some tensorial quantities on M || . In an 
adapted coordinates frame these may be easily constructed as follows 

9fiu = t M • t;, induced metric (first fundamental form ) 
a* = n* ■ d u tfj, second fundamental form 
jHj _ n « . q n j n ormal fundamental form 
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where = denote the tangent vectors to M associated with the chosen 
coordinate frame and the dot the standard scalar product in R n . 

The choice of an adapted coordinates frame is obviously not unique. An 
arbitrary coordinates transformation on M as well as a point dependent ro- 
tation of the normal vectors n l (x) transforms an adapted coordinates frame 
into an adapted coordinates frame. Whereas varying the choice of the coordi- 
nates x^ causes g^ v , a 1 and A l i to transform as tensors of M, the variation 
of normal vectors n l (x) by a rotation R kl (x) makes a 1 to transform as a 
SO (n — m) vector but Ajf as a SO(n — rn) gauge connection 

a% — ► n ik A kl n jl + ii ik dji jk . (4) 

The normal fundamental form A l j actually represents the connection induced 
by R n on the normal bundle of M, TM ± . 

The metric Gu, I,J— 1, ...,n of R n in the adapted coordinates frame 
y 1 } writes 

r _ (i,v + y k y l A kh A l u h y k Aj?\ 

Gu-[ ykA ik gu ) (5) 

where, introducing the matrix rj by 77^ = y l g tip a l pin the matrix 7 may be 
written as 

lm = g, p (l - 77)^(1 - r])t (6) 

The determinant \G\ of Gu coincides with that of the matrix 7, I7I, and the 
inverse of the metric tensor may be calculated as 



ryVpykj^ki fiij _|_ ykyl j^ik jS^jl^po 



" \ r.,vp„,k Aki Xij 1 „,k„,l Aik Ajl^.pa- j • \' ) 



3 The Perturbative Expansion 

We come now to the dynamical aspect of the problem. The quantization 
of the system is performed very easily and unambiguously in R n , before 
considering the constraint. In a cartesian coordinates frame {r 1 ; I — 1, n}, 
dynamics is described by 

n = ~d I d I + v c , (8) 



5 



where Vc is the potential realizing the constraint. Ti acts on wavefunctions 
if) G C 2 (R n ) normalized with the condition 

\)\ 2 dr n = I. (9) 

In an adapted coordinates frame {x^, y l ; /i = 1, m, % = m + 1, n} Hamil- 
tonian (||) takes the form 

n = -^^dj^iGi^dj + v c , (io) 

and the normalization condition (^) transforms to 

f \ip\ 2 \G\ ll2 dx m dy {n ' m) = 1. (11) 

Since we are looking for an effective dynamics on the submanifold M we find 
it natural to perform a similitude transformation in such a way that wave 
functions are correctly normalized in C 2 (M) instead of C 2 (R n ). The aim is 
achieved by 



1 1/4 

ji/T 



|G 


|l/4 


M 


1/4 



I 11/4 |/^|l/4 \ lz ) 

n_j \J I ni 

| G |l/4 n 

where |g| denotes the determinant of the metric g^ u induced on M. Con- 
sidering the explicit form (|7]) of the inverse metric G IJ and introducing 
dfj, = + iA^Lij/2, where Ly = —i(y l dj — y^di) are the SO(n — m) genera- 
tors, Hamiltonian fllPf ) takes the quite complicated form 



n 



d ^\ 1/2d *TZm ~ ouiiMuiiM ^br^^U + Vfc- (13) 



2 | 7 |l/4^l'l ^| 7 |l/4 2|(/| 1 /4| 7 |l/4 
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|l/4 



At this point, and only at this point, the constraint is imposed by con- 
sidering conditions CI and CI. Condition CI assures that the potential Vc 
may be replaced by its power expansion around the minimum y — 0. Con- 
dition C2 states that there exist an adapted coordinates frame in which Vc 
only depends on the normal coordinates y. Without loss of generality the 
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constant term of the expansion may be neglected and the quadratic term 
diagonalized by means of a point independent rotation in the normal space 

V c {y) = + a ijk y i y j y k + b ijkl y i y j y k y l + ... . (14) 

The scale of the proper frequencies id 1 has been readsorbed in the adimen- 
sional parameter e _1 . The smaller e the deeper is the minimum of Vc and 
the more the system is squeezed on the constraint surface. 

e appears as a natural perturbative parameter in the theory and, rescaled 
the normal coordinates by y — ► e 1 / 2 ?/, a perturbative theory may be set up 
by expanding Hamiltonian ( ]i~3|) in powers of e 

eH = + + e 3 /2#(3/2) + ^#(2) + + 

+e 5/2 a ijk y i y j y k + e 3 b ijkl y i y j y k y l + ... . 

The constants a^k, bij k i, ... appearing in the expansion of the confining 
potential are such that the second, the third and further terms of the right 
hand side of (|TJ]) are small compared to the first term and in this sense 
are e-dependent. In practical application they appear in the perturbative 
expansion as independent parameter so that, for example, e > f 2 aij k y l y^y k is 
not in general of order e 5//2 , its magnitude depending on the explicit form of 
the potential Vc- The zero and first order terms of expansion (|i~5D has been 
extensively discussed in Ref. [|J. 

In accordance with Heisenberg principle the zero order dynamics depends 
only on normal degrees of freedom, 

= 1 (_ a . a . + ^V 2 ), (16) 

describing a system of (n — m) uncoupled harmonic oscillators with frequen- 
cies cu m+1 , ...UJ n . 

More surprising results follow from the analysis of the first order term, 
H {1) = --^ [d, + iAj'Ly) g^g 1 ' 2 [d v + Z -A k jL kl ^j + Q«(x), (17) 

where the potential may be expressed in terms of the intrinsic scalar 
curvature R and the extrinsic mean curvature £ as 

Q^) = ±R{z)~?(x). (18) 
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Aside from the potential term Q^ l \ H^ 1 ' is proportional to the Laplace op- 
erator on M coupled to the motion in normal directions by means of the 
minimal interaction with the gauge field A^Lij/2. It is therefore resonable 
to expect that in a perturbative picture describes the effective dynamics 
on M. This has actually been found in Ref.[[|. The surprising result, un- 
expected and unrecoverable by means of a formal treatment of constraints, 
is that the effective dynamics is coupled with gauge fields and quantum po- 
tentials induced by the intrinsic and extrinsic geometrical properties of the 
constraint surface. The physical relevance of such a geometry-induced dy- 
namical structure has been recently discussed in Ref. 0, showing how this 
phenomenon is observable in the effective rotational motion of some simple 
polyatomic molecules. 

Since we are interested in a realistic description of a constrained micro- 
scopic system we never consider the limit e — > 0. e is a small but finite pa- 
rameter, its magnitude depending on the characteristics of the system under 
consideration. It is therefore very important to know the explicit expression 
of further terms of expansion (|15D in order to predict the spectrum of the 
system with an adequate precision. 

To evaluate the explicit expression of the generic term of the expansion 
( |ToD we start by observing that the first and second terms of Hamiltonian 
( |13|) may be rewritten solely in terms of 7^ and In |1 — 77 1 as 

1 m ia dih^^di , \.. = —-didi + 

2 Ly 1/4 ' ' ' Ly 1/4 2 

+\{didi In \l-rj\) + ^ m 11-771)0% In |1-77|), 
4 o 

and 

1 Irtl 1 / 4 1 
d ^ v h\ 1/2 d — = --V V^V + 

2M 1/4 N 1/4 " N 1/4 2 

+ \ (V M 7^V, In |1 - V \) + (V M In |1 - V \) (V.ln |l - V \) , 

where, denoting by V M the covariant derivative associated with the connec- 
tion induced on M, 

W, = W,+ % -A^L iy (19) 
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It is very convenient to introduce the matrices 

i$0 = (N + l)g^y k a% a ^...y^aX aN g^, (20) 

V(o) = 9^ '■ The expansion in e of r y^ v and In |1 — r/| may then be computed 

as 



7 



E e 

Af=0 



JV/2 \w 
'I(N) ' 



(21) 



In 1 1 — 77 1 



00 N/2 

E^r 



7V=1 



iV 



(22) 



The evaluation of the iV/2-order term of the perturbative expansion flT5|) 
reduces so to a matter of simple algebra yielding 



)(N/2) 



iV > 2, where the potentials Q^ 2 ) may be written as 

QW = ~tr [9^] tr [9^] - -tr [9^9^] 
Q( 3 / 2 ) = -tr [^77] tr [77^77] - -tr [178^] - -tr [V^V^ 



(23) 
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and for iV > 4 



Af-3 



Q 



1 



(AT/2) 



AT- 2 

E 

K=0 



-tr 



?7 9i?7 



tr 



N-K-2 



dill 



1 

tr 

4 



- V -tr 
^ 4 

AT=0 ^ 



N-if-3 



V.77 + (N — K — 3)^77 



77^ 2 dir]dir] 



+ 



+ 



+ E E 

ft-=0 L=l 



i7 £-1 (V„Tj) 



tr 



AT-A^L-3 



(V„fj) 



It is remarkable that the perturbative expansion of the Hamiltonian of 
the system is completely written in terms of the induced metric g pu , the sec- 
ond fundamental form a 1 , its first and second covariant derivative V p a l , 
V CT Vpa;* and the normal fundamental form A 1 ^. A l j appears in the pertur- 
bative expansion only by means of the minimal coupling ([19]). 
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4 Spectrum and Effective Dynamics 



To evaluate the spectrum of the system we proceed now by means of the 
standard Raleigh-Schrodinger perturbation theory. We identify with 
the unperturbed Hamiltonian and the rest of expansion (|T5|) with the per- 
turbation V e , 

eH = H {0) + V e . (24) 

As stated in the previous section, represents a system of (n — m) 
uncoupled harmonic oscillator. We denote by XAfiv) its eigenf unctions, hav- 
ing collected the harmonic oscillator quantum numbers n m ±i, n n in the 
multiindex Af = (n m+ i, ...,n n ). The corresponding eigenvalues are given by 
j uj l {rii+l/2). The spectrum is degenerate every time the frequencies 
uj % satisfy linear conditions in the integer field. The zero order eigenfunctions 
corresponding to an energy are given by 

^M > (x,y) = M%)xAf{y) (25) 

and present an infinite degeneracy given by the presence of the arbitrary 
function of x, (pj^(x), beside that labelled by the multiindex corresponding 
to the energy 

The first order correction to E^°\ E^\ is obtained by diagonalizing the 
perturbation on degenerate states, that is by solving the Schrodinger equation 

U m (j){x) = £ (1 V(z), (26) 

where the Hamiltonian TC E ° is obtained by bracketing the order e term of 
V e , between the harmonic oscillator states corresponding to -E^ ) and 

<f)(x) is a vector wavefunction having as component the 4>n(x) with energy 
E^. The explicit expression of H is 

H E(0) = - 2i7a(10„ + ^)^ 1/2 (1^ + *40 + Q {1 \x) + Q {1 \x) (27) 
with 

A.^-A^U^ (28) 

Q {1) = \sTA%J^ ((LijLki) ~ (LyHLu)) , (29) 
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where (Lij) and (LijL^i) denote the matrices obtained by bracketing Ljj and 
LijLki between the harmonic oscillator states corresponding to and 1 
is the identity matrix with the dimension of the degenerate space. 7i E{0) 
appears as a free Hamiltonian on the constraint surface coupled with the ge- 
ometry induced gauge fields (|28|) and the potentials ( |18D and ([^). Equation 
( |2"U| ) have therefore to be interpreted as the Schrodinger equation describing 
the effective dynamics induced on the constraint surface. Note that for a 
surface £ embedded in the three-dimensional euclidean space R 3 the Hamil- 
tonian (|2T| ) reduces to (|2]). We do not comment anymore on this fact referring 
to d, | for details. 

Denoting by /C the quantum numbers labelling the eigenfunctions of H. E(0) 
and supposed the degeneracy to be completely removed, the eigenvalues of 
eTL are evaluated by means of the standard formula 



This allows to calculate the spectrum of the system with an arbitrary accu- 
racy as a power series in the parameter e. 

5 Particle Constrained on a Circle 

As a very simple but nontrivial example we consider a particle constrained 
to move on a circle embedded in R 3 by an harmonic potential. This al- 
lows to illustrate some peculiarities of constrained quantomechanical sys- 
tems which are systematically ignored in formal treatments. Let therefore 
be c : [0,2irR] — > R 3 , c(x) = (Rcos(x/R),Rsm(x/R),0) the embedding 
map of the circle in the three-dimensional euclidean space R 3 . The curve 
is parametrized by the arclength x, so that its tangent, normal and bi- 
nomial may be immediately evaluated as t(x) = (— sm(x/R), cos(x/R), 0), 
n(x) = (cos(x/R),sm(x/R),0) and b(x) = (0,0,1). Every smooth assign- 
ment of an orthonormal basis of the normal space to c in x may be obtained 
by rotating normal and binormal by a point dependent angle w(x) 




E<® + (AT,/C|P e |AT,/C)+ 



|(A/-,/C|n|AT',/C')| 2 



+ ... .(30) 



n 2 = cosu> n + sinu> b, 
n 3 = — sin w n + cos w b, 



(31) 
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where w(0) = w(2ttR) + 2ttz, z being an integer. The induced metric, the 
second fundamental form and the normal fundamental form of the embedding 
read 



9n = 1, 

2 1 

a u = — cos w, 
1 

a?i = — — sin w. 



(32) 



'ii 

1 23 



R 

Af = -w. 

2 3 

7] is so written as rj = ^-cosw — ^-sinw, whereas the covariant derivative 
( |T§D on c reads V x = d x — iwL 2 z- The direct calculation shows that V x ?7 = 0, 
so that the whole perturbative expansion (^) may be easily evaluated as 



9 v z / 2 

1,„ 1 



H {1) = -~{d x -iwL 23 " 2 



2 v x « y gi?2 . 



2 „,3 



(33) 



Af-2 



ffW J ) = (JV - 1) [ y - cosw - y - sinw) 
\ R R I 



The spectrum of the system may now be calculated by means of perturbation 
theory. As in the general case the infinite degeneracy of the zero order 
states is removed by solving the Schrodinger equation (p6|) for the effective 
dynamics on the constraint surface. In correspondence to the zero order state 
labelled by the harmonic oscillator quantum numbers (^2^3) the effective 
Hamiltonian on the circle writes 

= -\ (d x - zw(L 23 )) 2 + \ ((L 2 23 ) - (L 23 ) 2 ) w 2 - JL (34) 

where angled brackets denote again expectation values between harmonic 
oscillator states corresponding to the energy = u) 2 (ri2 + 1/2) + uj^{n 3 + 
1/2). 

Let us now discuss the physical meaning of the function w(x). If the 
confining potential is symmetric, that is u 2 = u) 3 , it is possible to choose the 
harmonic oscillator basis is such a way that L23 is diagonal. The effective 
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potential ((L?, 3 ) ~~ (L23) 2 )w 2 /2 vanishes then identically and (£23)^ appears 
then as a pure gauge field in the theory and may be removed by a different 
choice of normal coordinates y 2 ,y 3 . The effective dynamics on the circle and 
the whole perturbative expansion result then considerably simplified. On 
the contrary, if the confining potential is not symmetric, u 2 ^ u 3 , (L 2 s) = 
0, (£23) 7^ and a different choice of normal coordinates would cause the 
confining potential to be x-dependent. The effects produced by w(x) may 
therefore not be eliminated. Pictorially we may assimilate our model to a 
particle moving in a ring with a small ellipsoidal section. The function w(x) 
describes then how the section wraps up when moving along the ring. If the 
ring's section reduces to a circle (uj 2 = uj 3 ) then does not matter how the 
wrapping is done and we are always reconduced to the case w — 0. On the 
contrary the wrapping produces observable effects when the ring's section is 
not circular (oj 2 7^ uj 3 ). 

The case uj 2 = uu 3 being straightforward we concentrate on uu 2 7^ uj 3 . The 
effective Hamiltonian describing the dynamics on the circle reduces then to 



1 1 

2 x + 2 




» 2 -^-(35) 



Different choices of the wrapping function w(x) produce a completely dif- 
ferent effective dynamics. An arbitrary positive, everywhere finite smooth 
potential may be reproduced by a suitable choice of w. 

The simpler case we may consider is that in which the potential wraps 
us uniformly, say z times, w(x) = zx/R. The effective Schrodinger equation 

■ k t 

on the circle is then immediately solved by <P(n 2 ,n 3 ),k(x) = e l R x /\f2itR, k any 
integer, and 




(36) 



Finer corrections to the spectrum may be evaluated by going over in pertur- 
bation theory. This matter not being particulary interesting for this model 
is postponed to the next example which is physically more significative. The 
remarkable fact we learn from this example is that the realization of the 
constraint, that is the particulary form of the confining potential Vc charac- 
terizes the spectrum and the effective dynamics of the constrained quantum 
system. Such informations are completely lost within a formal treatment of 
the constraint. 



13 



6 Particle Constrained on a Sphere 
(The Rigid Diatom) 

As a second example of a constrained quantomechanical system we consider 
the motion of a particle on a sphere embedded in R 3 . To get some physical 
intuition on what we are dealing with let us consider a diatom. Aside from 
effects connected to the geometric phase [lQl, the effective Hamiltonian de- 



scribing the rotovibrational degrees of freedom of the molecule is written in 
the adiabatic approximation as 

h 2 ( d 2 d 2 d 2 \ 

where r = (x, y, z) is the relative position of the nuclei, \i the reduced mass 
of the system and Vbo is the Born-Oppenheimer potential. It is usual to 
assume that 

~ Vbo presents a deep minimum in correspondence of the molecular equi- 
librium length r , 

- Vbo depends only on the relative distance of the nuclei r = |r| and not 
on the orientation of the molecule in space. 

That is, Vbo behaves as a potential confining the motion from the ro- 
tovibrational configuration space R 3 to the sphere of radius r . Hamiltonian 
(|3?D describes therefore a constrained quantomechanical system in the sense 
we specified before, erf. (§). 

In order to adapt coordinates we introduce the usual angles 9 and <fi 
parametrizing the sphere and the normal coordinate y = ^-y^(r — r ). 
/ = fiiQ is the momentum of inertia of the diatom and to the frequency 
introduced by the Born-Oppenheimer potential, on = y ^-^P |r=r - The 
adimesional scale factor e = Ti/Iuo appears naturally in the definition of y 
once the zero order energy hu is factorized from the Hamiltonian. From 
most diatoms e is a very small parameter, e ~ 10~ 2 — 10~ 4 , and, as our 
notation anticipates, gives a measure of the rigidity of the molecule. The 
metric of R 3 in the adapted coordinates frame reads 

1 /(l + e 1/2 ?/) 2 0' 

(l + ^ 2 y) 2 sm 2 6 | , (38) 

V o 



G 



i j 



hue 
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while the Born-Oppenheimer potential writes as 



V BO = tiu {^y 2 + ay 3 + by 4 + ••••) , (39) 

where a = ah 1 ^ 2 / fi 3 ^ 2 uj 5 ^ 2 and b = bh/fi 2 u 3 , a and b being the usual spectro- 
scopical parameters. The rigidity parameter h/Iu plays therefore the rule of 
the parameter e we introduced in section 3 for a generic constraint^. A com- 
parison of ( 38 ) with equations fl5|) and (|6]) allows to write down immediately 



the induced metric and the second fundamental form on the sphere 

1 fl 

huje\0 sin 2 



9^ = T-Al J^)' ( 4 °) 



V 2 /l 



and 



a ^ = -J^{o sin 2 e)- (41) 

The codimension of the constraint surface being one, the normal fundamental 
form vanishes identically. The embedding of the sphere in R 3 is standard and 
the perturbative expansion H/huo = H^ + eH^ + ... + ay 3 + by 4 + ... is easily 
evaluated as 

H^ = \{-d 2 y+ y% 

1 / 1 „ 1 



H {1) = -- (—d e sm9d e + -^d 2 ) 
2 \srn9 sm^ 9 J 

= (-l) N (N-l)y N - 2 H^\ 



(42) 



The zero order Hamiltonian takes into account the vibrational motion of 
the diatom. is the angular momentum operator describing the effective 
rotational dynamics as that of a spherical top. The rest of the perturbative 
expansion reproduces the Dunham expansion [IT| of the nonrigid rotator 
taking into account rotovibrational interactions. The rotovibrational spec- 
trum of the diatom finds therefore a very natural interpretation in terms of 
constrained quantum mechanics. Given up the classical idea of constraint 

1 Note that the normal coordinate y appears as already rescaled. 
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(e — > 0) the rotovibrational structure appears naturally as a consequence of 
the physical structure of the constraint. The algorithm we present in this pa- 
per gives an automatic way to compute rotovibrational interactions and may 
result usefull in the analysis of rigid polyatomic molecules ||. For the mo- 
ment we conclude by evaluating the spectrum of the particle constrained to 
the sphere. Introduced creation/destruction operators relative to the normal 
coordinate the computation results algebraic in nature and may be performed 
to an arbitrary order in perturbation theory by means of computer algebraic 
manipulation. We report here the expansion of the energy £ n j to the third 
order in perturbation theory 

+^\ 1 4i l ( l + !)] + 4 W + !)] (« + I)" ~ I W + I)]') + (43) 



+ ( -6 - 3e 3/2 a a 2 



15 .„\ / 1 



n 



2 



4 / V 2, 

Replacing the values of e, a and b Eq. (fl3|) reproduces the standard expression 



of rotovibrational energies of diatoms [12 



7 Concluding Remarks 

The reduction of the motion of a quantomechanical system from its configu- 
ration space to a submanifold is by no means unique, in the sense that it is 
impossible to perform this operation by completely disregarding the motion 
in the directions normal to the constraint surface. Quantum mechanics is a 
field theory, after all, and the wave function of the system keeps on explor- 
ing the whole configuration space even if squeezed on the constraint surface. 
When the system is in an eigenstate of the confining potential we can obtain 
an effective Hamiltonian describing the dynamics on the constraint surface. 
As is clearly illustrated by examples of sections 5 and 6 this effective dynam- 
ics, described by H Ei °\ depends both on the specific normal eigenstate and 
on the explicit form of the confining potential. In any case the eigenvalues of 
the effective Hamiltonian give only the first order corrections to the spectrum 
of the system at finite e. An accurate description requires also the analysis of 
the interaction between the motion normal and along the constraint surface. 
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The perturbative expansion flI5|) we present in this paper takes into account 
this effect. 

Performing the limit e — > after subtracting the divergent zero order 
energies of the system produces a well defined description of the motion on 
the constraint surface. Nevertheless we consider this operation as artificial, 
the physical nature of the constraint being in the small but finite value of e 
(cfr. the discussion of the diatom). The whole perturbative expansion (tOf) is 
therefore necessary to characterize the dynamics of the constrained system. 

Aside from its conceptual importance, the perturbative expansion (|i~5|) 
may reveal of practical importance in the analysis of electrons confined on 
arbitrary surfaces and wires as well as in the analysis of polyatomic molec- 
ular spectra. The effective rotational dynamics of some simple polyatomic 
molecules has already been considered in Ref . || demonstrating the physical 
relevance of the induced gauge structure and quantum potentials (|2lf), (|T8]) 
and (p9[). Our hope is that expansion flip] ) may serve as an unifying tool in 
understanding the fine structure spectra of rigid polyatomic molecules. 
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